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We demonstrate that large class of PT-symmetric complex potentials, which can have isospectral 
real partner potentials, possess two different superpotentials. In the parameter domain, where the 
superpotential is unique, the spectrum is real and shape-invariant, leading to translational shift 
in a suitable parameter by real units. The case of two different superpotentials, leading to same 
potential, yields broken PT-symmetry, the energy spectra in the two phases being separated by a 
bifurcation. Interestingly, these two superpotentials generate the two disjoint sectors of the Hilbert 
space. In the broken case, shape invariance produces complex parametric shifts. 
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I. INTRODUCTION 



The fact that two different superpotentials can give rise to same potential in supersymmetric quantum mechanics 
(SUSY-QM) is well-known The superpotential is deformed by an additive function, which leads to the construc- 
tion of isospectral Hamiltonians having a number of physical applications. Complex potentials, particularly the ones 
having parity (P) and time reversal (T) invariance have attracted significant attention in the recent literature 0. 
Appearance of real eigenvalues in certain parameter domain for these complex PT-symmetric potentials, have led to 
considerable interest in this class of not so well-studied spectral problems. It has been found that the spectrum is real, 
O" 1 when the wavefunction respects PT-symmetry, whereas complex eigenvalues, paired by complex conjugation are real- 
ized when the wavefunctions do not respect the above symmetry. As a function of certain potential parameter(s), the 
. energy eigenvalues show bifurcation, when the spectra transits from real to complex values. Detailed discussions and 
' the progress in this field can be found in the following references Q and . A number of models have been studied, 
both numerically and analytically for illustrating the above structure in [5|, and subsequently utilised in It was 
observed that complex PT-symmetric potentials, under suitable parametrization, can be isospectral to real potentials 
0. A generalized Poschl- Teller- type complex potential, of the above form, has been constructed by Ahmed 3, which 
reveals all the above mentioned features explicitly. Interestingly, it was found that SUSY only yields the real spectra, 
where wave-functions respect PT-symmetry. Explicit solutions of differential equation led to the complex branch of 
the spectra under different parametric conditions, where the solutions do not respect PT-symmetry. It is then natu- 
ral to ask, what is the precise relationship between SUSY and PT-symmetry and if the aforementioned result is generic. 



In this paper, it is shown that the presence of complex parameters in the superpotential can lead to the real- 
ization of the PT-symmetric complex potential through two different superpotentials. As compared to the case of 
j_j • isospectral deformation, in the present case the superpotentials do not differ by an additive function. The difference 
in superpotentials appear through the complex constant parameters. In the broken PT phase, two superpotentials 
are present leading to the same potential, whereas in the unbroken phase the correspondence between the potential 
and superpotential is unique. This is possible only for complex potentials and does not manifest in the case of the real 
ones. The two different parameter domains explain the observed bifurcation separating the unbroken PT case from 
the broken one. In the latter phase, shape-invariance leads to complex translational shifts in the relevant parameters, 
whereas in the former case the translation is in real units [9[ . Very interestingly, in case of complex eigenvalues, half 
of the Hilbert space is generated by one superpotential, whereas the other one generates the remaining part of the 
wavwfunctions . 

The paper is organized as follows. In the following, we provide a brief introduction to SUSY-QM and pro- 
ceed to construct the potentials from a general superpotential, in section III, starting from Ahmed's potential. 
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After illustrating the nature of PT-symmetry in the broken sector, we derive the spectra, in both broken and 
unbroken phases, through shape-invariance. We then list examples of a large class of potentials exhibiting the above 
characteristics. In section IV, a number of solvable complex potentials are constructed through the above procedure, 
which are not PT-symmetric. These shape-invariant potentials constitute complexified Poschl- Teller-type fl9| and 
Coulomb-type potentials, which undergo complex translation in the appropriate parameter domain. Finally, we 
conclude after pointing out a number of interesting directions in which the present investigation can be advanced. 



II. SUPERSYMMETRIC QUANTUM MECHANICS 

Supersymmetric quantum mechanics Q], fiol. ITU fl2T| interrelates the spectra of two different Hamiltonians, H±, which 
can be written in the factorized form [13[, where H_{x) = A^A and H+{x) — AA^ (h = 1 = 2m): 

H ± (x) = -^ + V ± (x). 



In terms of the superpotential W(x), the potentials are 

dx 



V ± = W*(x)±^; (1) 



here. 



"'-(-s+ww). 



and W(x) 



ipo (x) being the ground-state cigcnfunction of iJ_ (x) . 



1 dip (x) _ 
ipo(x) dx 



It can be straightforwardly shown that, if tp~(x) is an eigenfunction of Atp~(x) is an eigenfunction of 

H + with the same eigenvalue, except for the ground state of H_ defined as Aip (x) = 0. Similarly if (x) is an 
eigenfunction of H+, then A^ip~(x) is an eigenfunction of Taking into account the unpaired ground state ipo(x), 
the respective eigenvalues are related as, 

The corresponding eigenfunctions satisfy, 

#t"(aO = [E- +1 }^A^- +1 , 
and ip- +l {x) = [E+]-iA^~. 

It was shown by Gendenshtein that if two isospectral potentials satisfy the relation 

V+(x;a ) = V-(x;a 1 )+R(a 1 ), (2) 

where do is a set of parameters for a given pair V±, a\ = /(do); and R(ai) is independent of x, then one can construct 
a hierarchy of Hamiltonians, 

o2 s 
H s = - +V-(x;a s )+J2 R ( a k) 
x fc=i 



dx 2 



V+(x;a s _i) + y^R(a k ) 



k=l 
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with ground state energy 



s 



£ s = ]Ti?(a fe ). 



k=i 



On identifying that H 1 = H + and H° = H-, it is found that the spectrum of if_ is given as 




fe=i 



Isospectral potentials satisfying Eq.(2) are called shape-invariant [HI]. Using the relation ip n (x]a ) oc 
A'(x; ao)ipn( x 'i a i) it can be shown that 



ip n (x; a ) oc A^(x; a )A*(x; a{) A^(x; a„-i)ip (x; a n ). 



Therefore for shape-invariant potentials, one can solve the eigenvalue problem purely algebraically and obtain the 
spectrum and corresponding eigenfunctions. Detailed analysis of the properties of PT-symmetric quantum mechanical 
systems exhibiting supersymmetry has been carried out by Mostafazadeh [l5[ and Plyushchay et al. [l6j], including 
the construction of appropriate norm [l5j . The question of norm in general pseudo-Hermitian Hamiltonian systems 
has been answered recently [17j . 



Complex PT-symmetric potentials are known to have real eigenvalues, despite not being Hermitian in the usual 
sense. Further, this realness of spectra survives only over a specific range of parameter values, beyond which 
PT-symmetry is broken. In the broken PT range one observes complex conjugate spectra, as a function of certain 
potential parameters, with the corresponding states connected by PT-operation. If such a potential is shape- invariant, 
SUSY can be applied to algebrically obtain the spectra and corresponding eigenfunctions. 

Here, we start with a general superpotential, which leads to both real and complex spectra. Under cetain 
parametrization, the superpotential is unique, leading to real eigenvalues and corresponding wave- functions. For 
a different parametrization, two different superpotentials yield the same potential, leading to complex conjugate 
eigenvalues with corresponding eigenfunctions related by PT-operation. As a result, we arrive at the supersym- 
metric parameter condition for phase-transition of the spectrum from real to complex conjugate values owing to 
spontaneous breaking of PT-symmetry. The experimental observation of phase-transition in a PT-symmetric system 
has recently been reported by Guo et al. [l8l |. 

To start with, we consider Ahmed's potential [8j], which is given as, 



Explicit solution of the Schrodinger equation shows that it has both real and complex-conjugate spectra, depending 
upon the relation among the real parameters V\ and V2, given as |Va| < V\ + \ and | V2 1 > V\ + j, respectively 

In order to see the role of supersymmetry and PT-symmetry, we start with the complex superpotentials 



where A, B, C PT are real constant parameters. Instead of repeating C PT in both the coefficients above, in general, 
one can start with different parameters. We have chosen them to be same to arrive at Ahmed's potential. The general 
case will be discussed afterwards. The V-(x) corresponding to Eq.(5) is, 



III. CONSTRUCTION OF COMPLEX PT-SYMMETRIC POTENTIALS 



V(x) 



V\sech 2 (ax) — iV2sech(ax)tanh(ax) . 



(4) 



W^ T = (A± iC PT )tanh(ax) + {±C PT + iB)sech(ax), 



(5) 




[(A ± iC PT )(A ± iC PT + a)- (±C PT + iBf] sech 2 (ax) 
-i(±iC PT - B) [2(A ± iC PT ) + a] sech(ax)tanh(ax), 



(6) 
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which in general may not be PT symmetric. To be PT-symmetric, the coefficient of the first term must be real and 
that of the second term must be purely imaginary. Implementing these conditions, one arrives at a unique relation: 



C P1 [2{A-B) + a] =0. 



(7) 



This is the most important result of the present paper, which we now study in more detail. 



Condition. 1: If C PT = 0, then 

Wpt(%) = W rea i(x) = Atanh(ax) + iBsech(ax), 



(8) 



which gives 



V ± (x) = V-(x) = - (A(A + a) + B 2 ) sech 2 (ax) + i (B(2A + a)) sech(ax)tanh(ax). 



(9) 



This matches with Ahmed's potential, with the identifications, 



Vi = A{A + a) + B 2 , 
and V 2 = -B(2A + a), 



and satisfy the reality condition |V^| < Vi + \. From Eqs. (1),(2), (3) and (9) we get the corresponding known 
spectrum 



E n = -{na-A) 2 , 



(10) 



modulo a constant term. One can also obtain the corresponding eigenfunctions: 

B 



ip n {x) oc {sech(ax)) a exp[— i — tan (sinh(ax))]P n 



A B 1 A _ 



[isinh(ax)] 



(11) 



using Eq.(4) as per the prescription in 0], which matches with those given in Q for real spectrum. 



Condition. 2: For C ^ 0, from Eq.(7), one gets A = B — ^, which when substituted in Eqs. (5) and (6) yields 



Wp T (x) = W±{x) = (A±iC PT )tanh(ax)+ ±C PT + i (a + ^) sech(ax) 



(12) 



and 



V ± (x) = Vl{x) 



2A(A + a) -2(C PT ) 2 



a 



sech 2 (ax) 



2A(A + a) + 2(C 



PT\2 



sech(ax)tanh(ax) . 



(13) 



From the above we get 



Vt = 2A{A + a) ~2(C P1 Y + 



a 



4 ' 



and V% = — 



2A(A + a) + 2(C pl ) 2 + — 
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which satisfiy the condition | | > V% + |- It needs to be emphasized that both the superpotentials in Eq.(5) leads 
to the same PT-symmetric complex potential, corresponding to broken PT symmetry. We can obtain their spectrum 
through shape-invariance, 

B± = 2n (A ± iC PT ) a + (na) 2 , (14) 



which has complex conjugate pairs and each member of a complex conjugate pair can be connected to one of the two 
superpotentials in Eq.(12). Same can be said about the corresponding eigenfunctions, 



tpn( x ) ^ (sech(ax)) 



±(A±iC PT ) 



exp 



C 



PT 



tan 1 (sinh(ax)) 



,24 + * 



isinh(ax)] , 



(15) 



which match with the previous results Q, under SUSY parametrization. They are related to each-other through PT 
transformation. 

As mentioned earlier, one can start with two different parameters, say C PT and D PT . Then we arrive at 
two conditions, 



D PT = 



2BC 



PT 



2A + a' 

and [{2A + a) 2 - (2B) 2 ] C PT = 0, 



for the potential to be PT-symmetric. Then if C PT ± 0, we have ±(2A + a) = 2B and D PT = ±C PT , which leads 
to the same results for the broken PT case. If C PT = 0, then A and B are unrestricted, giving the same results as in 
the unbroken PT case. 



The key observation about the potentials obeying the above property is that the odd parity part of the po- 
tential must have a purely imaginary coefficient and the even part must have a purely real one. With these 
constraints, one can construct a large class of PT-symmetric potentials, which are shape-invariant, starting from the 
known real shape-invariant potentials. Below, in table I, we list the superpotentials, the condition for which spectral 
bifurcation takes place leading from a PT-symmetric ground state to one which does not respect the same. 
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Table I: List of shape-invariant PT-symmetric potentials, with their respective superpotentials and parametric 
conditions for spectral bifurcation. The corresponding energies and ground state wave- functions are also shown. 
(Here we have taken C PT = C for simplicity.) 



W{x;A,B,C) 


Condition 


Cases 




Energy(#„) 


Ground-state 


(A ± iC) 
xtanh(ax) 


C(2A + a) 
= 


C = 


— A(A + a)sec 2 (ax) 
+2iBtanh(ax) 


A 2 - (A - no) 2 

B 2 i B 2 


— B 

(sech(ax)) " exp(-i^x), 
A > 




7? ' (A-na) 2 


A B 
~ rL {A±iC) 






(A 2 + C 2 ) sec 2 (ax) 
+2iBtanh(ax) 


(A±iC) 2 
-(A±iC -no) 2 

B 2 
{A±iC) 2 
B 2 

' (A±iC-nu) 2 


(sech(ax) ) ~ ( - A±lC " > 
xexp (—gr^i {iA± C)s) 
A > 


-(A±iC) 
xcoth(ax) 


C(2A - a) 
= 


C = 


A(A — a)csch 2 {ax) 
—2iBcoth{ax) 


A 2 - (A + na) 2 

B 2 B 2 


(sinh(ax))~ exp {— i^rx) , 
~^<A<0 




A 2 > (A+na) 2 


+i (A±iC) 




A=f 


- (A 2 + C 2 ) csch 2 {ax) 
—2iBcoth(ax) 


(A±lC) 2 

-(A±iC + na) 2 

B 2 
(A±iC) 2 
B 2 

' (A±iC+na) 2 


(sinn(ax)) a 
xexp f-^afga {iA ± C)x\ , 


(A ± iC) 
xtan(ax) 


C{2A - a) 
= 


C = 


A(A — a)sec 2 (ax) 
+2iBtan(ax) 


-A 2 + (A + na) 2 

B 2 i B 2 


(cos(ax))~ exp {— ijrx) 




A 1 (A+na) 2 




"' B 

^ l {A±iC) 






-(A 2 +C 2 ) sec 2 (ax) 
+2iBtan(ax) 


-{A±iC) 2 
+(A ±iC + na) 2 

B 2 
(A±iC) 2 

+ (A±iC+na) 2 


(cos^ax))" 
xexp (^~ A2 + C 2 (iA±C')x S j 


(A ± iC) 
xcot(ax) 


C(2A + a) 
= 


c = o 


A(A + a)csc 2 (ax) 
+2iBcot(ax) 


-A 2 + (A - na) 2 

B 2 i B 2 


(sin(ax))~~ exp {— i^x) , 




IF " r {A-na) 2 


+i B 




^ = -1 


-{A 2 + C 2 )csc 2 (ax) 
+2iBcot(ax) 


-(A±iC) 2 
+ {A±iC -na) 2 

B 2 
(A±iC) 2 
B 2 

^ (A±iC-na) 2 


(sin(ax)) 
xexp (- A2 ^ c a {iA ± C)xj , 
< A < 



The complex conjugate spectra appear due to spontaneous breaking of PT-symmetry, the specific parametrization 
condition C PT ^ can be identified as the SUSY criteria for broken PT. This is different from the analytic parameter 
criterion, where PT is unbroken over a range of parameters. In contrast, here we have a discreet condition C PT — 
for the same. Furthermore, broken PT results in the bifurcation of the corresponding Hilbert space in terms of 
two distinct superpotentials; the corresponding eigcnfunctions map into each-other under PT operation. They 
correspond to a unique complex potential VI (x), which is still PT-symmetric as required. Apart from a, the Ahmed 
potential has two independent parameters owing to the condition A = B — ^ , just like the corresponding potential 
for unbroken PT following C PT = 0. Same can be seen in case of all the potentials listed above. As mentioned 
earlier, the non-uniqueness of the superpotential is well known from isospectral deformation, which incorporates an 
additional function obeying Bernoulli's equation Here, the non- uniqueness arises parametrically. 

IV. NON-PT-SYMMETRIC POTENTIALS 



In the previous section, we have seen how a generic isospectral complex potential, under certain parametric 
constraints, can demonstrate a link between PT-symmetry and SUSY-QM. We now study shape-invariant complex 
potentials which need not have PT-symmetry. Since these potentials are shape-invariant, one can obtain their spectra 
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analytically. 

We start with a minimal complexification of the well-known Poschl- Teller fl~oj ] potential: 

U(x) — U a sech 2 (ax) + Ubcsch 2 (ax), (16) 



by considering two superpotentials 



Wi(x) = Atanh{ax) + iBcoth(ax), 
and W2{x) = iAtanh(ax) + Bcoth(ax). (17) 



The respective isospectral potentials are, 



U^(x; A, B) = -A(A =F a)sech 2 (ax) - B(B ± ia)csch 2 (ax), 
and C/^(:r;A, 5) = A(A±ia)sech 2 (ax) + B(B a)csch 2 (ax). (18) 



Both of them are shape-invariant since, 



U+{x; A, B) = [/f (x; A - a, B + ia) + (A + iB) 2 -{A + iB - 2) 2 , 
and U£(x;A,B) = U~(x; A + ia,B - a) + (iA + B) 2 - (iA + B - 2) 2 . (19) 



This immediately yields their respective spectra: 

E„ = 4na(A - na) + iABna, 
and E 2 = Ana(B - na) + iAAna. (20) 



The separations between two successive energy levels AE n = E n — E n -i are given respectively as 

A E\ = -4a(A - 2na - a) + i2Ba, 
and A E 2 = —Aa(B - 2na - a) + ilAa. (21) 



The corresponding eigenfunctions are, 



ip n ( x ) ^ [cosh(ax)] a [sinh(ax)] a P n " 2 ' a 2 (sinh(axj) , 
and ip n (x) oc [cos/i(air)] a [szn/i(aa;)] Q P„ ° 2 ' ° 2 (sinh(axj) . (22) 



It is straightforward to check that, when the above potential is real, one finds the corresponding energies to 
be E^ eal — 4na(A + B — na), with energy spacings AE^ eal = — 4a(A + B — 2na — a). The corresponding 

-A -B. I^B.-! -A-l\ 

eigenfunctions are (a;) oc [cosh(ax)] a [sinh(ax)] a P n a 2 ' " 2 (sinh(ax)). 

Due to the minimal complexification the original real spectra turns into a complex spectra with the real part 
mimicking the original spectra. Interestingly, the imaginary part of the energy displays equispaced harmonic 
oscillator-like spectra. The non-uniqueness of the original spectra arising from infinite possible combinations of 
parameters A and B, giving the same (A+B), is also removed. As in both the examples of complexification, only the 
parameter that got multiplied with i appears in the equispaced imaginary spectra, the complexification process does 
not differentiate between the two functions appearing in the superpotential. 

From Eqs.(22), we see that the powers appearing in the wave-functions become imaginary. In case of the 
first wave- function in Eqs.(22), this completely removes the singularities of the function sinh(ax) for the whole range 
of B. Similarly, the singularity of the function cosh(ax) is removed in the second case for the whole range of A. 



8 



Thus wave-functions for these complexified potentials are normalizable over a greater range of parameters than those 
belonging to their real counterpart [20j] . 



One can construct a complex Coulomb type potential starting with the superpotential, 



which results in the potentials 



1 CV 

W{r) = — +0, 
r 



a(a ± i) i2al3 ~ 



(23) 



If the coefficient of K af3 is constant(say e), then one can compare V- with the radial Coulomb potential 

1(1 + 1) e 
Vc(r) = —3 -, 



where I and e carry the conventional meanings. The shape-invariance condition, 

" 1 I 



V + (r;a) = V- (r; a + i) + <y 2 



a 2 (a + iy 



is satisfied and one finds, 



a 2 (a + inf- 



(24) 



which is similer to the Hydrogen spectra. Starting from the superpotential in Eq.(23), the complex eigenfunctions 
can be obtained 11( . For example the ground-state wavefunction is given by 



^o(r) ocr ia exp- /3r , (25) 
which is well-behaved over the whole range of the real parameter a, unlike the well known real counterpart. 



V. CONCLUSIONS 



In conclusion, for a large class of potentials, it is found that a given PT-symmetric complex potential can be 
realized from two parametrically different superpotentials. In the parameter domain, where the superpotential is 
unique, the Hamiltonians yield real eigenvalues and the latter leads to the broken PT phase. In this case, the 
two superpotentials yield two disjoint parts of the Hilbert space of the Hamiltonian and shape-invariance leads to 
complex shifts in the appropriate potential parameters. Higher order supersymmetry has been applied to construct 
new classes of PT-symmetric potentials [21], which may also reveal the features found here. Our procedure also 
yields non-PT-symmetric potentials with complex eigenvalues, which can be obtained through shape-invariance. 
This feature may manifest in many-body systems of the Calogero-Sutherland type [13, HH, [24|. This dynamical 
system is currently under investigation. Finally, we have studied here translation-type shape-invariance. It is worth 
investigating, through our method complex supersymmetric potentials where shape invariance arises through scaling 
of the parameters. 

Acknowledgment: We acknowledge useful discussions with Prof. R. Dutt, and are thankful to Profs. A. 
Mostafazadeh and Plyushchay for bringing to our notice a number of related references. 
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